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ABSTRACT 

We report on the analysis of hard-state power spectral density function (PSD) of GX 339-4 
down to the soft X-ray band, where the disc significantly contributes to the total emission. At 
any luminosity probed, the disc in the hard state is intrinsically more variable than in the soft 
state. However, the fast decrease of disc variability as a function of luminosity, combined with 
the increase of disc intensity, causes a net drop of fractional variability at high luminosities 
and low energies, which reminds the well-known behaviour of disc-dominated energy bands 
in the soft state. The peak-frequency of the high-frequency Lorentzian (likely corresponding 
to the high-frequency break seen in active galactic nuclei, AGN) scales with luminosity, but 
we do not find evidence for a linear scaling. In addition, we observe that this characteristic 
frequency is energy-dependent. We find that the normalization of the PSD at the peak of 
the high-frequency Lorentzian decreases with luminosity at all energies, though in the soft 
band this trend is steeper. Together with the frequency shift, this yields quasi-constant high 
frequency (5-20 Hz) fractional rms at high energies, with less than 10 percent scatter. This 
reinforces previous claims suggesting that the high frequency PSD solely scales with BH 
mass. On the other hand, this constancy breaks down in the soft band (where the scatter 
increases to ~ 30 percent). This is a consequence of the additional contribution from the 
disc component, and resembles the behaviour of optical variability in AGN. 
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1 INTRODUCTION 

Black hole X-ray binaries (BHXRB) spend most of their time in 
a quiescent state, which is occasionally interrupted by outbursts of 
activity. As the outburst evolves the source passes through differ¬ 
ent accretion states and both the variability and spectral properties 
change (e.g. Miyamoto et al. 1992; Belloni et al. 2005; Remillard 
& McClintock 2006). The variability properties of these sources 
have been extensively studied in the past, revealing the existence 
of trends as a function of accretion state and/or luminosity (e.g. 
Zdziarski et al. 2004; Munoz-Darias, Motta, & Belloni 2011; Bel¬ 
loni, Motta, & Munoz-Darias 2011). However, due to the sensitiv¬ 
ity limits of the detectors employed (particularly the Proportional 
Counter Array, PCA, on the Rossi X-ray Timing Explorer, RXTE), 
these studies have been mostly limited to energies > 3 keV. At these 
energies the broad-frequency-band (e.g. in the typical frequency 
range ~ 0.1 - 60 Hz) aperiodic variability progressively decreases 
as the outburst evolves, with the fractional root-mean-square (rms) 
variability amplitude (e.g. Nandra et al. 1997; Vaughan et al 2003; 


Ponti et al. 2004) dropping from values of tens of percent in the 
canonical hard state, down to a few percent in the canonical soft 
state (Belloni et al. 2005; Munoz-Darias, Motta & Belloni 2011). 
During the hard state the fractional rms shows an “inverted spec¬ 
trum” (decreasing with energy), which switches to a “hard spec¬ 
trum” (increasing with energy, at least above ~5 keV) during in¬ 
termediate and soft states (e.g. Zdziarski et al. 2004; Belloni et al. 
2011, as illustrated in the inset of Fig.|4j. These behaviours are usu¬ 
ally ascribed to variations of the relative contribution of the main 
X-ray spectral components (e.g. Gierlihski & Zdziarski 2005). In¬ 
deed, over the past years, the general idea was that the strong broad- 
frequency-band variability characterizing the hard state was associ¬ 
ated with the Comptonized hard X-ray emission (power law) from 
the corona (or the base of a radio-emitting jet, e.g. Zdziarski et al. 
1998; Markoff, Nowak, & Wilms 2005; Droulans et al. 2010), since 
this component dominates the X-ray spectrum during this state. On 
the other hand, the disc-thermal component, known to be variable 
on the time scales of the outburst, was thought to be intrinsically 
constant on shorter (than about hours-minutes) time scales (e.g. 
Gierlihski & Zdziarski 2005). Given that the disc dominates the 
spectrum up to energies of ~ 10 keV during the soft state, it would 
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then be responsible for the drop of broad-frequency-band variabil¬ 
ity amplitude, and for the change of shape of the fractional rms 
spectrum. 

The way these variability properties extrapolate to soft X-ray 
energies (< 3 keV) during the hard state is currently poorly 
known. Since, during the hard state the disc component gives non- 
negligible contribution in the soft X-ray energy band ( i 1 keV; 
Tomsick et al. 2008; Kolehmainen, Done, & Diaz-Trigo 2014), 
to some extent we would expect to observe behaviours similar to 
those characterizing the disc during the soft state. However, the re¬ 
cent use of data from the EPIC-pn detector onboard XMM-Newton 
(sensitive down to ~0.3 keV), for the spectral-timing analysis of 
BHXRBs, unveiled the existence of high levels of intrinsic (on time 
scales of i 1 s) disc variability in the hard state (Wilkinson & Ut- 
tley, 2009). Though this result has been obtained studying a single 
dataset for each of the two sources GX 339-4 and Swift J1753.5- 
0127, it might represent a general property of BHXRBs. Thus, ex¬ 
tending timing studies down to soft X-ray energies is crucial in or¬ 
der to gain a complete understanding of the disc-corona variability 
processes. 

Another important property of X-ray variability in BHXRBs is that 
the characteristic time scales vary during the outburst. During the 
hard state all the variability components (commonly modeled with 
Lorentzians; Belloni, Psaltis, & van der Klis 2002) in the power 
spectral density function (PSD) shift towards higher frequencies 
as the luminosity increases, while keeping tight correlations be¬ 
tween each other (e.g. Belloni et al. 2005; Done & Gierlihski 
2005). This indicates strong dependence of X-ray variability on 
the luminosity/mass-accretion rate of the source. The shift of the 
characteristic frequencies as a function of mass-accretion rate is 
a general property of accreting systems. Indeed, as in BHXRBs, 
a scaling (consistent with being linear; McHardy et al. 2006; Ko- 
erding et al. 2007) also links the “high-frequency break” (Edelson 
& Nandra 1999) and the mass-accretion rate in active galactic nu¬ 
clei (AGN). However, despite all these variations in both character¬ 
istic frequencies and broad-frequency-band variability amplitude, 
the high-frequency PSD (i.e. above the peak-frequency of the high- 
frequency Lorentzian) remains remarkably constant during the en¬ 
tire outburst (Gierlihski, Nikolajuk & Czerny, 2008), indicating 
little or no-dependence on the luminosity and the accretion state 
of the source. In analogy with BHXRBs, also in AGN the high- 
frequency PSD (above the high-frequency break) does not show 
any clear dependence on the luminosity/mass-accretion rate (e.g. 
O’Neil et al. 2005; Zhou et al. 2010; Ponti et al. 2012). 

In this paper we aim at investigating the evolution of the PSD as 
a function of luminosity/mass-accretion rate, extending the analy¬ 
sis down to the soft X-ray energy band. Eor this reason, we study 
the hard-state PSD of GX 339-4, the best monitored BHXRBs. We 
used data collected with XMM-Newton, and, in order to obtain a 
broad band coverage, we considered also simultaneous RXTE ob¬ 
servations. The data selection and reduction are presented in Sect. 

while the PSD estimate and analysis procedures are described in 
Sect. (3] Given the different, frequency-dependent behaviours out¬ 
lined before, which characterize the PSD at energies i 3 keV, 
we will first study the broad-frequency-band variability proper¬ 
ties (Sect. 14. H4.2b . and then more specifically focus on the high- 
frequency part (Sect. l43t . Results are discussed in Sect.j^ 



Hardness 

Figure 1. The hardness-intensity diagram of GX 339-4 as obtained from 
all the archived RXTE observations relative to the 2004 (light grey dots) 
and 2009 (dark grey triangles) outbursts. The luminosity is computed in the 
range 3-30 keV, the hardness is computed as the ratio between the 13-30 
keV and 3-13 keV flux. Colored points mark the position of the analyzed 
XMM-Newton and simultaneous RXTE observations in the diagram. The 
labels refer to the nomenclature used in the paper to indicate the 2009 (low 
luminosity), 2004 (medium luminosity), and 2010 (high luminosity) obser¬ 
vations. The error bar on the left side of the plot reports the uncertainty on 
the HID normalization due to the uncertainties on the distance and the BH 
mass of GX 339-4. 

2 DATA REDUCTION 

GX 339-4 has been observed 18 times by XMM-Newton (as of May 
2014). In this paper we focus on all the archived observations with 
exposures >10 ks, which caught the source in a hard state, i.e four 
observations in total. In order to further increase the spectral band 
coverage we considered also all the available simultaneous RXTE 
observations. Details about the analysed data sets are reported in 
Tabled 

Among the four selected XMM-Newton observations, two 
(16-18 March 2004) were carried out during consecutive satellite 
revolutions. Hereafter, we will show results obtained by combining 
these two observations together, and we will refer to them as one 
single observation. The remaining two observations belong to the 
same outburst of the source, but they have been carried out one year 
apart (in 2009 and 2010). 

We used EPIC-pn (Striider et al. 2001) and Proportional Counter 
Array (PCA; Jahoda et al. 2006) data, which ensure high effective 
area over a broad spectral range. 

The XMM-Newton data reduction was done using the XMM Sci¬ 
ence Analysis System (SAS vl3.5), and the latest calibration files 
(CCF, as of May 2014). The analyzed data sets are all in Tim¬ 
ing mode. We followed standard reduction procedures, and applied 
Rate Dependent PHA (RDPHA) and X-ray loading (XRL) correc¬ 
tions (Guainazzi & Smith 2013; Guainazzi 2014, and references 
therein). We selected time intervals free from background proton 
flares (GTI, the final effeclive exposures are listed in Table [T]l and 
considered only evenfs wifh PATTERN < 4. The source counts 
were extracted from the region which includes 90 percent of the 
total collected counts. For all the observations this region falls be¬ 
tween RAWX columns 31 and 45. 

Being all the observations in Timing mode, pile-up is expected to 
be an issue only when the count rate exceeds ~800 cts/s (in the 0.7- 
10 keV range, Guainazzi et al. 2014). This threshold is well above 
the average count rate registered during the 2004 and 2009 obser¬ 
vations considered here. However, the 2010 observation exceeds 
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Table 1. Log of the analyzed XMM-Newton and simultaneous RXTE observations: (1) observation ID; (2) date; (3) detector used for the analysis and 
corresponding observing mode (for the RXTE observations the Proportional Counter Units - PCU - coadded in the analysis are also reported); (4) effective 
exposure (the XMM-Newton effective exposures refer to the duration of the observation after removal of high soft proton background, the percentage of this 
exposure affected by telemetry drop outs is reported within parentheses) (5) mean count rates in the energy bands used for the analysis (i.e. bands 0.5-1.5 keV 
- soft - and 2-9 keV - hard — for XMM-Newton, and band 10-30 keV - very hard - for RXTE). 


XMM-Newton 

(1) 

ObsID 

(2) 

Date 

(3) 

Instmment/Mode 

(4) 

Exp (drop outs) 

[s] 

(5) 

Mean count rate (soft/hard) 
[cts/s] 

0204730201 

2004-03-16 

Epic-PN/Timing 

104200 (20.7 percent) 

90/108 

0204730301 

2004-03-18 

Epic-PN/Timing 

71000 (21 percent) 

98/117 

0605610201 

2009-03-26 

Epic-PN/Timing 

32220 (0.2 percent) 

43/61 

0654130401 

2010-03-28 

Epic-PN/Timing 

33630 (23.7 percent) 

357/434 

RXTE 

ObsID 

Date 

Instmment/Mode 

Exp 

Mean count rate (very hard) 




Is] 

[cts/s/PCU] 

90118-01-05-00 

2004-03-16 

PCU 0,2,3/Good Xenon 

2324 

55 

90118-01-06-00 

2004-03-17 

PCU 0,2/Good Xenon 

1660 

58 

90118-01-07-00 

2004-03-18 

PCU 0,2/Good Xenon 

6505 

60 

94405-01-03-00 

2009-03-26 

PCU 0,2/Good Xenon 

3055 

38 

94405-01-03-01 

2009-03-26 

PCU 2/Good Xenon 

3305 

37 

94405-01-03-02 

2009-03-26 

PCU 1,2/Good Xenon 

1807 

39 

95409-01-12-01 

2010-03-28 

PCU 2/Generic Event 

19195 

180 


this threshold. In addition, variability may cause pile-up during 
short time intervals even when the average count rate is below the 
threshold. Therefore we verified the pile-up level during each ob¬ 
servation using the SAS tool “epatplot”. In the 2010 and 2004 ob¬ 
servations, the observed single and double events fractions slightly 
deviate from the expected trend for data not affected by pile-up (re¬ 
spectively by ~1.2 and 12 percent in the 2-10 keV band for both 
observations), but this deviation is systematic at all energies. As 
noticed in Done & Diaz-Trigo (2010), these systematic deviations 
might be either due to the uncertainties related to the modeling of 
the singles and doubles fraction for the Timing mode, or might in¬ 
dicate that the data are actually slightly affected by pile-up. Thus, 
we tested our results against pile-up by repeating the analysis (de¬ 
scribed in Sect.O of the 2010 observation (which has the highest 
EPIC-pn count rate) on the event file obtained by excising the cen¬ 
tral three pixels (i.e. RAWX= 37 - 39), which reduces the excess of 
double events due to pile-up by a factor ~5 in the 2-10 keV band. 
We found results perfectly consistent with those presented in Sect. 
[4] the only difference being the expected reduction of variability 
power signal-to-noise (S/N) ratio in the PSD. 

For the RXTE data analysis we followed standard procedure|3 us¬ 
ing the HEASOFT software v6.13. For the timing analysis, we used 
either Good Xenon or Generic Event configuration modes, depend¬ 
ing on the availability (see Table [T](. To maximize the S/N ratio 
we extracted data from all the layers of the Proportional Counter 
Units (PCU) that were simultaneously and continuously switched 
on during each observation. On the other hand, the luminosities are 
computed from data in Standard2 configuration, considering only 
the top layer of PCU2. 

The position of the analyzed observations within the hardness- 


* The reduction procedures are described in the RXTE data reduction cook¬ 
book: http://heasarc.nasa.gov/docs/xte/recipes/cook_book.html 


intensity diagram (HID) is displayed in Fig.[T] This plot includes all 
the RXTE data relative to the 2004 (light grey dots) and 2009 (dark 
grey triangles) outbursts. Fluxes are computed in the 3-30 keV en¬ 
ergy range, assuming a simple power law plus a gaussian at the en¬ 
ergy of the Fe Ko- line. The hardness ratio is computed as the ratio 
between the 3-13 keV and 13-30 keV fluxes. The fluxes are con¬ 
verted into Eddington-scaled luminosities. Eddington luminosities 
are obtained by considering the currently available estimates and 
uncertainties on the distance and the BH mass of the source (i.e. 
~6-10 kpc, and ~6-10 Mq, respectively, Hynes et al. 2003, 2004, 
Munoz-Darias et al. 2008) and assuming a reasonable mean value, 
i.e. 8 kpc and 8 Mq. The error bar displayed on the left side of 
the plot highlights the uncertainty on the plot normalization, as due 
to the uncertainties on the distance and the BH mass. Note that, 
when considering the relative increase of luminosity, the relevant 
uncertainty is the error on the flux measurements, which however 
is very small (i.e. of the order of ~ 0.1 - 0.3 percent), thus it has 
not been reported in the plot. The analyzed observations show lu¬ 
minosities Li-iOkeVILedd ~0.007 (2009 observations), ~0.02 (2004 
observations), and ~0.07 (2010 observations). In the following we 
will refer to these observations respectively as “low”, “medium”, 
and “high luminosity”. 


3 DATA ANALYSIS 

For each observation we extracted light curves with time bins of 
1ms (we verihed that this choice of the time bin does not cause 
any spurious feature within the analysed frequency window, which 
might affect the analysis and results presented in this paper) in the 
three energy bands: 0.5-1.5 keV, 2-9 keV (using XMM-Newton 
data), and 10-30 keV (using RXTE data). Hereafter we will refer 
to these bands respectively as soft, hard, and very hard. Accord¬ 
ing to previous analyses of the same data sets (Kolehmainen et al. 
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Table 2. Results obtained from the fit of the PSDs of GX 339-4 with a 
model comprising three broad Lorentzians. Best-fit values of the peak- 
frequencies (in units of Hz) of the lowest, intermediate, and highest fre¬ 
quency Lorentzian (respectively, v/, y,, and v/,) are reported, in the soft (0.5- 
1.5 keV, S), hard (2-9 keV, H), and very hard (10-30 keV, VH) band, and at 
different luminosities. Errors are reported at 90 percent confidence level. 


Band 

Low L 

MedL 

V/ 

High L 

s 

0.0096 ± 0.0014 

0.042±0.002 

0.146±0.013 

H 

0.0092 ± 0.0014 

0.043± 0.003 

0.172±0.026 

VH 

0.0080 ± 0.0035 

0.043±0.004 

0.179±0.016 

V/ 

s 

0.18±0.04 

0.56±0.05 

0.166±0.008 

H 

0.22±0.02 

0.76±0.05 

0.161±0.008 

VH 

0.17±0.06 

0.55±0.08 

0.161±0.012 

V/, 

s 

1.15±0.15 

1.18±0.10 

1.85±0.15 

H 

1.88±0.15 

2.79±0.29 

2.89±0.11 

VH 

1.90±0.42 

2.58±0.44 

3.05±0.28 

T^/dof 

S 

80.6/42 

161.9/66 

21.4/25 

H 

95.9/42 

87.1/66 

60.3/25 

VH 

44.4/42 

89.4/30 

36.8/25 


2014; Plant et al. 2015), the chosen soft band is where the accretion 
disc gives the most significant contribution. In the other two energy 
bands the dominant component is the primary hard X-ray power 
law (this is shown also in Sect. ld.l.TT l. 

To study the evolution of the variability characteristics of GX 339- 
4 at different energies, we estimated the PSD within each of these 
energy bands. To this aim, the light curves were sampled over in¬ 
tervals shorter than the total exposure and of equal length within 
each observation. The length of the segments for the XMM-Newton 
light curves was chosen so as to exclude the gaps due to teleme¬ 
try drop-outs (which occur mostly when the target source is very 
bright, i.e. during the high and medium luminosity observations), 
while ensuring a sufficiently broad frequency coverage. For consis¬ 
tency, the RXTE light curves have been sampled according to the 
corresponding simultaneous XMM-Newton light curves (i.e. 21 s, 
59 s, and 500 s for the high, medium and low luminosity obser¬ 
vations, respectively). We estimated the PSD (e.g. Nowak et al. 
1999b) of each light curve segment, using routines implemented 
with IDL v6.4, and averaged them to obtain an estimate of the in¬ 
trinsic PSD for every observation and energy band. The level of 
variability power due to Poisson noise was determined fitting the 
high-frequency part (above 100 Hz) of the PSD with a constant. 
Then, the average PSDs have been Poisson noise-subtracted, and 
renormalized adopting the squared fractional rms normalization 
(Miyamoto et al 1991). We corrected the RXTE data for instru¬ 
mental dead-time (which produces an over-estimate of the intrin¬ 
sic Poisson noise level) using the dead-time model of Zhang et al. 
(1995, 1996). Since the variability power S/N ratio drops at high 
frequencies, we limited the analysis to frequencies i 50 Hz. 

The hard-state PSDs of BHXRBs are usually very complex. They 
contain several structures which can be modeled with multiple 
Lorentzian components (e.g. Nowak 2000, Belloni et al. 2002, 
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Figure 2. The PSDs of GX 339-4 in the three energy bands 0.5-1.5 keV 
(soft), 2-9 keV (hard), and 10-30 keV (very hard), for each observation 
at different luminosity. The arrows mark the hest-fit peak frequency of the 
low-frequency (v/) and the high-frequency (v*) Lorentzians for each energy 
band. 


Pottschmidt et al 2003). We followed this same approach for the 
modeling of the PSDs of GX 339-4. 


4 RESULTS 

The hard-state PSDs of GX 339-4 are shown in Figs. |2] and In 
Fig.E] we show for each observation, the PSDs at different ener¬ 
gies, while in Fig. we show for each energy band, the PSDs at 
different luminosities. 

All the PSDs show a similar shape, with two broad humps, whose 
intensity and characteristic frequencies clearly change as a func¬ 
tion of both luminosity and energy. These two broad humps are a 
characteristic feature of BHXRBs in the hard state (e.g. Nowak, 
Wilms, & Dove 1999a; Wilkinson & Uttley 2009; Bock et al 2011; 
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Figure 3. Overplotted are the PSDs of GX 339-4 at different luminosities, 
for each of the energy bands soft (upper panel), hard (middle panel), and 
very hard (lower panel). 


Grinberg et al. 2014; Cassatella et al. 2012). At frequencies lower 
than those covered here the variability power per unit logarithmic 
interval (vPy) is expected to drop in the hard state (band-limited 
noise). The chosen length of light curve segments allows us to ob¬ 
tain a good sampling of the broad-band noise of GX 339-4, down 
to the onset of this drop. We verified the absence of additional low- 
frequency broad structures by inspecting the PSDs obtained using 
longer segments of RXTE data. This confirmed that the analysed 
PSDs are band-limited, at least in the hard and very hard bands. 
We note that the soft band PSD shows the same structures as those 
seen at higher energies, therefore we assume that also the soft band 
PSD is band-limited. 

We used two Lorentzian components to model the two humps and 
more rigorously quantify their variations. A third Lorentzian is re¬ 


quired in all the fits to account for residual structure^. The best- 
fit centroid frequency and full width at half maximum have been 
used to estimate the peak-frequency (i.e. the frequency of the max¬ 
imum power per logarithmic interval; Belloni et al. 2002) of each 
Lorentzian (see Table [3 using the formula reported in Sect. 14.21 
The quality of the data is high enough, that our simple model, 
though providing a good description of each PSD, does not al¬ 
ways return sufficiently small values of the reduced (see Ta¬ 
ble [2j- We verified that this is mainly due to the presence of small 
residuals which are left unmodeled. These residuals are usually ac¬ 
counted for by increasing the number of Lorentzian components 
(e.g. Nowak 2000, Belloni et al 2005). However, comparing our 
results with those found in the literature, both for GX 339-4 (e.g. 
Nowak at al. 1999a, Nowak 2000) and other BHXRBs (e.g. Cas¬ 
satella et al. 2012, Grinberg et al. 2014) we conclude that three 
variability components are sufficient to obtain a good description 
of the hard-state broad-band noise of GX 339-4 within the anal¬ 
ysed frequency window (i.e. up to frequencies of 20 Hz). Our best- 
fit models are plotted in Fig.[2(solid lines) together with the single 
components (dotted curves). To avoid confusion, only the compo¬ 
nents associated with the very hard band PSD are shown in the 
plots. These best-fit PSD models were first used to study the vari¬ 
ations of the broad-frequency-band (between 0.002-20 Hz) frac¬ 
tional rms variability amplitude (Sect. iTTb . Then we used them to 
constrain the characteristic frequencies and study their variations 
as a function of energy and luminosity (Sect. lAa . Finally, we fo¬ 
cused on the variability properties of the high-frequency (between 
5-20 Hz) portion of the PSD (Sect. 14. 3t . 

4.1 Trends of broad-frequency-band fractional rms as a 
function of energy and luminosity 

We integrated the best-fit PSD models for each luminosity and 
energy band, over the frequency interval 0.002-20 Hz (broad- 
frequency-band fractional rnis), respectively for the high, medium, 
and low luminosity observations. Note that there is no evidence 
of additional low-frequency broad noise structures in the RXTE 
data (Sect. 0, and therefore these best-fit models describe the data 
well down to 0.002 Hz, at least in the hard and very hard bands. 
Given that the soft band PSD shows similar structures as observed 
at higher energies, we assume that the extrapolation of our best- 
fit model describes the data well down to very low frequencies 
also in this band. We also estimated the fractional rms associated 
with the low- and high-frequency Lorentzian (obtained integrat¬ 
ing the single components of the best-fit PSD models within the 
same broad-frequency-band interval). Results are reported in Table 
(3] The broad-frequency-band fractional rms, independently of lu¬ 
minosity and energy band, is always high ( i 30 percent), covering 
values typical of the hard state (Munoz-Darias et al 2011). In Fig. 
|4]we plot the broad-frequency-band fractional rms as a function 
of energy (within the three energy bands adopted throughout this 
paper). Note that the fractional rms of each Lorentzian shows the 
same energy-dependence as the broad-frequency-band fractional 
rms. For comparison we also report (inset of Fig. 0 the spectral 
shape of the 0.1-32 Hz fractional rms of GX 339-4 during a (high 
luminosity) hard and the soft state, at the energies covered by the 

^ The third Lorentzian fits the narrow peak of the low-frequency broad¬ 
band noise component of the high luminosity observation. In the other two 
observations, this Lorentzian fits an additional broad component at interme¬ 
diate frequencies (see Table[2). 
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Table 3. The fractional rms (e.g. Vaughan et al 2003) as obtained integrating 
the PSD best-fit models of each observation at different luminosities, and in 
each of the three bands 0.5-1.5 keV (soft), 2-9 keV (hard), 10-30 keV (very 
hard). The broad-frequency-band fractional rms and the fractional rms as¬ 
sociated with each Lorentzian are obtained integrating over the frequency 
interval 0.002-20 Hz. Finally we report the high-frequency fractional rms 
oblained integrating over Ihe interval 5-20 Hz. Errors are computed at the 
90 percent confidence level. 



Soft 

Hard 

Very Hard 

Broad-frequency-band 

HighL 

MedL 

Low L 

0.318±0.022 

0.492±0.034 

0.670±0.036 

0.376±0.024 

0.496±0.031 

0.600±0.036 

0.318±0.020 

0.374±0.024 

0.433±0.024 

Low-frequency Lorentzian 

HighL 

MedL 

Low L 

0.231 ±0.024 

0.363 ±0.037 
0.426 ±0.046 

0.251±0.026 

0.349±0.035 
0.384 ±0.041 

0.203±0.022 
0.261 ±0.027 
0.281 ±0.023 

High-frequency Lorentzian 

HighL 

MedL 

Low L 

0.187±0.020 
0.312±0.032 
0.386 ±0.030 

0.223±0.022 

0.240±0.028 

0.335±0.031 

0.189±0.019 

0.223±0.021 

0.245±0.026 

5-20 Hz 

HighL 

MedL 

Low L 

0.075±0.007 
0.109±0.010 
0.132± 0.010 

0.118±0.011 
0.138± 0.011 
0.147± 0.012 

0.103±0.010 

0.112±0.011 

0.108±0.010 


RXTE PCA (adapted from Belloni et al. 2011). 

Our main results can be summarized as follows: 

• Above ~2 keV, where the power law component dominates the 
source emission, we always observe a decrease of broad-frequency- 
band fractional rms as a function of energy. Therefore such a de¬ 
crease of variability with energy appears to be an intrinsic property 
of the power law component. This might be produced by variations 
of the spectral index, with a high-energy pivot point. Extrapolat¬ 
ing this trend onto the soft energy band, we see that the broad- 
frequency-band fractional rms in this energy band is lower than 
expected. This is consistent with the presence of a disc-component, 
less variable than the power law. 

• Comparing the decrease of broad-frequency-band fractional 
rms within each energy band as a function of luminosity, we find 
that the largest drop is observed in the soft energy band (a factor 
~2. 1, as compared to ~ 1.6 and 1.4 in the hard and very hard energy 
bands). In Sect. l4. l.Tl we show that, though the intrinsic fractional 
rms of the disc-component is high in the hard state, the fraction 
of variable disc flux in the soft band decreases fast (faster than the 
power law component), thus causing the observed drop of soft band 
variability as the luminosity increases. 

• In both the medium and low luminosity observations, the 
broad-frequency-band fractional rms spectrum has an overall “in¬ 
verted” shape (decreasing with energy), typical of the hard state 
(e.g. Zdziarski et al. 2004; Belloni et al. 2011; inset of Fig. |4ll. 
However, during the high luminosity observation, the large drop 
of broad-frequency-band fractional ms observed in the soft band 
causes a shift of the peak of variability power towards intermediate 
energies (i.e. in the hard band). This appears to be the onset of a 



Figure 4. Variations of the broad-frequency-band fractional rms spectrum 
with luminosity. Different symbols refer to the energy bands: soft (black 
dots), hard (red triangles), and very hard (blue squares). The inset shows the 
typical spectral shape of the 0.1-32 Hz fractional rms of GX 339-4 during 
a (high luminosity) hard and soft state observation in the RXTE bandpass 
(adapted from Belloni et al. 2011). 



Figure 5. The broad-frequency-band fractional rms as a function of disc- 
to-total (absorbed) flux ratio in the soft band. 

transition towards a “hard spectrum”, typical of intermediate and 
soft states (inset of Fig. 0. 

In order to characterize the nature of the drop of broad- 
frequency-band fractional rms observed in the soft energy band at 
high luminosity, and investigate how this can be associated with an 
increasing contribution of disc-thermal emission in this band, we 
carried out fits of the XMM-Newton energy and covariance spec¬ 
tra (e.g. Wilkinson & Uttley 2009; Uttley et al. 2011; Uttley et al. 
2014). 

4.1.1 Is the disc variable or constant during the hard state? 

Previous analyses of the same (or part of the) data sets analysed in 
this paper (Wilkinson & Uttley 2009; Kolehmainen et al. 2014) re¬ 
vealed the presence of a significant contribution from the accretion 
disc-thermal emission in the soft band. Thus we used the XMM- 
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Newton data to measure the disc-thermal emission-to-total flux ra¬ 
tio in the soft band and verify whether the drop of fractional rms is 
related to spectral variations in this band. We fit the XMM-Newton 
spectra between 0.5-10 keV using a simple model for the contin¬ 
uum (i.e. tbabs*[diskbb+nthcomp], with the Comptonization high 
energy cut-off fixed at 100 keV and the seed photon temperature 
tied to the disc temperature), and excluding the range of energies 
4-7 keV, dominated by the Fe Ko- line. We also excluded the 
energies around the edges of the response matrix, since residuals 
at these energies might be an artifact of uncorrected XRL effects 
and/or charge transfer inefficiency (Kolehmainen ef al. 2014). The 
confribufion from the disc emission strongly depends on the cold 
absorption column density parameter, Nh- To be more conserva¬ 
tive, we let this parameter free. Hence, all the uncertainties stated 
below account for the uncertainty on N/j. Nevertheless, the derived 
best-fit values range between Nh ~ 5 - 6 x 10^* cm“^, in agree¬ 
ment with the expected values (Dickey & Lockman 1990; Kong 
et al. 2000). In Fig. we plot the broad-frequency-band fractional 
rms as a function of disc thermal emission-to-total (absorbed) flux 
ratio in the soft band, for the different luminosities. A net increase 
of the disc relative contribution to the total flux as the luminosity 
increases is observed. This increase coincides with the net decrease 
of fractional rms, and the relation between the two is consistent 
with being linear (the best fit slope being -1.26 ± 0.24). Thus the 
observed drop of broad-frequency-band fractional rms at soft X- 
ray energies is strongly linked to the progressive increase of a disc 
component, in analogy with what would be expected from a com¬ 
parison with disc-dominated states. 

In order to produce a drop of the fractional rms, the disc com¬ 
ponent should be less variable than the power law. Thus we esti¬ 
mated the fractional rms intrinsic to the disc and the power law 
components. To this aim we computed the broad-frequency-band 
covariance spectra (with respect to the 0.5-10 keV reference band) 
of the XMM-Newton data (see Wilkinson & Uttley 2009; Uttley 
et al. 2011; Uttley et al. 2014), which allow us to determine the 
fraction of variable disc emission. Note that, strictly speaking, the 
covariance spectrum gives the spectral shape of the components 
whose variability is linearly correlated with the variability in the 
reference band. However, for intrinsic coherence close to unity (as 
in this case, see also De Marco et al. 2015), the covariance spec¬ 
trum is equivalent, though with smaller error bars (Wilkinson et 
al. 2009; Uttley et al. 2014), to the rms spectrum, which gives the 
spectral shape of all the variable components. We fit both the co- 
variance and the energy spectra with the same model used above, 
and with the Nh parameter tied to the value obtained from the fit 
of the energy spectra. Fig. shows the ratios to the best-fit power 
law, once the normalization of the disc-component is set to zero. 
All the covariance spectra show an excess at soft X-ray energies. 
This demonstrates the presence of a variable disc component at all 
the luminosities spanned by the analysed hard-state observations. 
The intensity of the excess in the covariance spectra is lower than 
in the energy spectra in the medium and high-luminosity obser¬ 
vations, indicating that some fraction of the total disc emission is 
constant or that the disc is less variable. We used the best-fit mod¬ 
els to derive estimates of the intrinsic fractional rms of the disc and 
the power law components (for each of the two components this is 
computed as the ratio between the absorbed flux in the covariance 
and in the energy spectrum). These are reported in Table |4] In ad¬ 
dition, we report the fractions of variable/constant power law/disc 
flux contributing to the total soft band flux. 

We observe that the disc is intrinsically variable at all lumi¬ 
nosities. Its fractional rms is comparable with that intrinsic to the 


Table 4. The fractional rms of the disc (in the soft band) and the power law 
(in both the soft and the hard band) components, and the fractions of the 
total soft band flux due to variable/constant disc/power law emission (Dyar, 
Pvar and Dconsu Pconst)- These estimates have been obtained by fitting the 
broad-frequency-band (0.002-20 Hz, 0.017-20 Hz, and 0.04-20 Hz, respec¬ 
tively for the low, medium, and high luminosity observations) covariance 
spectra of GX 339-4 and its XMM-Newton energy spectra with the model 
described in Sect. i4m All these estimates have been derived considering 
the absorbed fluxes. Errors are reported at 90 percent confidence level. 



Low L 

MedL 

HighL 

Disc [0.5-I.5keV] 

fractional rms 

0.65±0.20 

0.40±0.03 

0.20±0.03 

^var 

0.17±0.03 

0.15±0.01 

0.09±0.01 

^ const 

0.09±0.08 

0.24±0.01 

0.36±0.02 

Power law [0.5-1.5 keV] 

fractional rms 

0.68±0.02 

0.47±0.01 

0.34±0.01 

Pvar 

0.51 ±0.04 

0.29±0.01 

0.19±0.01 

P const 

0.23±0.03 

0.32±0.01 

0.37±0.01 

Power law [2-9 keV] 

fractional rms 

0.59±0.02 

0.42±0.01 

0.33±0.01 


power law component (at low luminosities it is consistent with be¬ 
ing as variable as the power law), and much higher than typically 
observed in soft states (a few percents; Munoz-Darias et al. 2011). 
However, the disc fractional rms decreases with luminosity faster 
than the power law fractional rms. Combined with the variation of 
the relative contribution from the disc and the power law to the to¬ 
tal soft band flux (Fig. O, the result is a net increase (by a factor 
of ~4) of constant disc emission in the soft band (Dconst) as a func¬ 
tion of luminosity. In comparison, the constant power law emission 
in the soft band (Pconst) increases only by a factor ~1.6. Thus we 
conclude that the drop of fractional rms in the soft band is driven 
by the net increase of constant disc flux in this band. It is worth 
stressing that all these considerations refer to the broad-frequency- 
band, while a more detailed treatment of the high-frequency (i.e. 
above the peak-frequency of the high-frequency Lorentzian) PSD 
properties is deferred to Sect. 14.3 1 

4.2 Testing the correlation between the characteristic 
frequencies of the PSD in the soft band 

Fig.lH reports the best-fit peak frequency of the low-frequency 
Lorentzian, v/, plotted against the high-frequency one, v* (these two 
frequencies are marked by arrows in the plots of Fig.|2]l. The peak 
frequency of each Lorentzian is computed using the formula (e.g. 
Belloni et al. 2002): 

V= ^l^+A2 

where vg is the centroid frequency and A is the half-width-at-half- 
maximum (HWHM) of the Lorentzian. Note that flat-top, broad¬ 
band noise components are best-fit by Lorentzians with vq con¬ 
sistent with zero. This is the case for most of our fits. However, 
in some cases, the high-frequency Lorentzian fit required a more 
peaked profile, with vo different from zero. The error on the peak 
frequency is computed from the 90 percent confidence contours of 
the Vo and A parameters. Both v; and V;, show the characteristic 
shift towards high frequencies as a function of luminosity. While v; 
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Figure 6. Data-to-model ratios of the energy (red squares) and covariance spectra of the low {left panel), medium {middle panel), and high luminosity {right 
panel) observations. The model used to fit the data is tbabs * [diskbb + nthcomp], and the normalization of the diskbb component has been set to zero to make 
the plots. All the covariance spectra have been computed using the 0.5-10 keV range as the reference band. The broad-frequency-band covariance spectra 
(black circles) are obtained integrating over the entire frequency range covered by each data set up to frequencies of 20 Hz (i.e. 0.002-20 Hz, 0.017-20 Hz, 
and 0.04-20 Hz respectively for the low, medium, and high luminosity observations). The high-frequency covariance spectra (light grey triangles) are obtained 
integrating over the frequency range 5-20 Hz. 





Figure 7. The best-fit peak frequencies of the high- and low-frequency 
Lorentzians. The colors refer to the energy bands: soft (black), hard (red), 
and very hard (blue). The grey points are the frequencies of the high- and 
low-frequency Lorentzian components as obtained by Belloni et al. (2005) 
from the analysis of all the hard-state RXTE observations of GX 339-4 
during the 2002/2003 outburst. The dashed curve is a linear best-fit to the 
Belloni et al. (2005) data points, while the dotted curve is the same linear 
model (with only the intercept left free to vary) fit to the soft band data 
points derived in this paper. 

shifts by a factor of ~20, V;, shifts by a factor of ~1.6, for a varia¬ 
tion in luminosity of a factor ~ 10. The black dots show that this 
same shift is observed also in the soft energy band. In addition we 
observe that v;, is systematically and significantly (at > Scr confi¬ 
dence level) offset towards lower frequencies in the soft band. Note 
that this offset is driven by a variation of A, indicating that the high- 
frequency Lorentzian is systematically narrower in the soft band, 


than at higher energies. It is also worth noting that, if the high- 
frequency PSDs (e.g. 1.5 times above v;,) are fitted with a simple 
power law, a significant steepening of the slope is observed in the 
soft band of the high-luminosity observation. This resembles the 
behaviour observed in AGN in the disc-dominated optical band 
(e.g. Mushotzky et al. 2011). On the other hand, v/ is consistent 
with remaining constant among the different energy bands of the 
same observation. Indeed, for each luminosity, the low-frequency 
Lorentzians in the three energy bands are co-aligned in the plots 
of Fig. (2] while an offset of v* in the soft band is clearly observed 
(see also Fig.[8l panel a, and Sect. l4.3l for a description of this fig¬ 
ure). Note that the peak frequency of the additional third Lorentzian 
does not show any significant energy-dependence (see also Table 
IS. These results suggest that the high-frequency PSD behaves in a 
different manner in the soft band than at higher energies. This as¬ 
pect will be explored in more detail in Sect. 14.31 
In order to compare our results with those reported in the literature, 
we overplot in Fig.|3(grey points) the characteristic frequencies as¬ 
sociated with the high- and low-frequency Lorentzians (in the en¬ 
ergy band £ = 3.8 - 15.3 keV) as obtained by Belloni et al. (2005) 
from the analysis of all the RXTE observations of the 2002/2003 
outburst of GX 339-4 (only the values relative to the hard state of 
the source are shown in the plot). The two main Lorentzians of the 
hard- and very hard-bunA PSDs analysed in this paper (red and blue 
points in Fig.|7]( follow the same trend and span the same range of 
characteristic frequencies as those of the 2002/2003 RXTE data. 
This trend seems to characterize also the Lorentzians in the soft- 
band PSD, despite the significant offset. To show that the correla¬ 
tion between v/ and v;, is preserved in the soft band, we fit the Bel¬ 
loni et al. (2005) data points with a linear model (in the log-linear 
plot of Fig.|3the model is represented by the dashed curve). Then 
we fixed the slope to the best-fit value, and letting the intercept free 
to vary, we fit the soft band data points (dotted curve in Fig. |7](. 
Though the number of data points is small, we obtain a very good 
description of the v/ vs. v;, trend in the soft band, thus in agreement 
with the idea that the two Lorentzians are linked by the same linear 
relation as observed at high energies. According to these fits, the 
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Figure 8. The scaling as a function of LjLEdd (where L is the 3-30 keV 
luminosity) of: (a) the peak-frequency of the high-frequency Lorentzian, 
va; (b) the normalization of the PSD at the peak-frequency of the high- 
frequency Lorentzian, ; (c) the high-frequency (5-20 Hz) fractional rms. 
The dashed lines are the best-fit models. In the plots, the colors refer to the 
energy bands: soft (black), hard (red), very hard (blue). 

offset of V;, as measured in the soft band and at higher energies is of 
a factor ~ 2. 


4.3 The high-frequency PSD of GX 339-4 

We now focus on the high-frequency PSD of GX 339-4 (between 
5-20 Hz, i.e. above the peak of the high-frequency Lorentzian). 
In the analysed datasets we observe that the peak of the high- 
frequency Lorentzian varies with luminosity by a factor of ~L6 
(Sect. O, as also shown in Fig. [8] (panel a), where we plot v/, 


as a function of (note that L is computed in the 3-30 keV 

range). A fit with a linear model in the log-log space returns a best- 
fit slope of 0.26 ± 0.11, 0.17 ± 0.06, and 0.19 ± 0.16, respectively 
in the soft, hard, and very hard band. We note that these slopes are 
consistent with being the same, but consistently and significantly 
different from a slope of 1. However the soft band has a signifi¬ 
cantly lower intercept (differing by a factor ~ 0.7 ± 0.2), this being 
a consequence of the energy-dependence of v;, (Sect. l4^ . 

We also looked at the variations of the normalization of the 
PSD at the peak-frequency of the high-frequency Lorentzian (W,,, 
in units of Hz^9 as a function of luminosity. We derived and 
its related uncertainties (at the 90 percent confidence level) from 
the best-fit PSD models. Results are shown in Fig. [8] (panel b). A 
fit with a linear model in the log-log space gives best-fit slopes 
of -0.81 ± 0.13, -0.52 ± 0.08, and -0.42 ± 0.20, respectively in 
the soft, hard, and very hard band (dashed lines in Fig. panel 
b). These results point to the existence of an inverse-power scaling 
relation between Ay^ and Once again we note a consistent 

behaviour between the hard, and very hard bands, and a significant 
deviation in the soft band. Indeed, in the soft band this trend appears 
steeper than in the other bands. 

Finally, we computed the fractional rms in the 5-20 
Hz frequency range. Results are reported in Table [3] and plot¬ 
ted in Fig. [8] (panel c). Despite the significant variations of the 
high-frequency Lorentzian as a function of luminosity, the high- 
frequency fractional rms (above the peak-frequency of the high- 
frequency Lorentzian) in the hard and very hard bands is consistent 
with being remarkably constant. This is due to the combined effect 
of the frequency shift and the decrease of normalization of the high- 
frequency PSD with luminosity. These combined variations pro¬ 
duce a constant high-frequency fractional rms. Indeed, fitting the 
data with a constant we measure a scatter around this best-fit value 
of ~ 10 and ~ 4 percent (respectively in the hard and very hard 
band), i.e. of the order of the corresponding average uncertainty 
(Fig.m middle and lower panel). On the other hand, the constancy 
of the high-frequency fractional rms with luminosity breaks down 
in the soft energy band (Fig.O upper panel). If fitted with a con¬ 
stant model, the increases to ~ 30 percent. Instead, a fit with a linear 
model in the log-log space yields a best-fit slope of -0.24 ± 0.08, 
and the scatter around the model is of ~ 3 percent. This is due to 
the fact that the normalization of the high-frequency PSD decreases 
with luminosity faster in the soft band than at higher energies. 


5 DISCUSSION 

In this paper we investigated the luminosity-dependence of the 
characteristic frequencies and variability power of GX 339-4 in the 
hard state. These properties have been extensively studied in the 
past, through the use of RXTE data, which cover energies > 3 keV. 
Our work aims at extending this analysis down to the soft X-ray 
energy band (£ > 0.5 keV), where the contribution of the disc 
thermal emission becomes non-negligible. For this reason we used 
XMM-Newton and simultaneous RXTE observations, covering the 
energy band 0.5-30 keV. The selected observations belong to two 
different outbursts of the source. Nonetheless, the variability prop¬ 
erties of BHXRBs are extremely similar from outburst to outburst 
(Munoz-Darias et al. 2011, Grinberg et al. 2014), thus justifying 
our approach of considering observations taken during different ac¬ 
tivity periods. 
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5.1 Disc variability and the shape of the 

broad-frequency-band fractional mis spectra 

The drop of fractional rms that characterizes the transition from the 
hard to the soft state, is ascribed to the emergence of a disc com¬ 
ponent, which, in the soft state, is more intense and less variahle 
(with typical fractional rms values of few percents) than the power 
law. The fractional rms spectrum is observed to change accordingly, 
from an “inverted” (in the hard state) to a “hard” (in the soft state) 
spectral shape. 

Extending our analysis of the hard-state variability properties 
of GX 339-4 down to the soft X-ray band, we observe a drop of 
fractional rms that correlates with the progressive increase of disc 
flux relative to the total flux in that band. This resembles the be¬ 
haviour seen in soft states, although to a much lower extent (indeed 
the disc fractional rms is much higher and the drop less intense 
than in soft states, see inset of Fig. |4l(, and is expected if the disc 
has lower fractional rms than the power law. This raises the ques¬ 
tion as to whether or not the disc in the hard state has similar timing 
properties as in the soft state. Analysing the energy and covariance 
spectra of GX 339-4, we observe intrinsic broad-frequency-band 
disc variability at all the luminosities spanned by the hard-state ob¬ 
servations. The presence of a variable disc in a hard state obser¬ 
vation (corresponding to our medium luminosity observation) was 
first reported by Wilkinson & Uttley (2009). We detect high levels 
of disc variability (~ 20 - 65 percent) in the three observations, 
comparable to the variability intrinsic to the power law component, 
and significantly higher than observed in soft states. However, the 
disc fractional rms decreases with luminosity faster than the power 
law fractional rms. At the same time, the fraction of soft band flux 
due to the disc (power law) component increases (decreases) with 
luminosity (Fig.O- The net effect is a strong increase (a factor of 
~4) of constant disc-flux fraction in the soft band, which causes the 
observed drop of fractional rms. It is worth noting that, without the 
contribution from the disc component, the fractional rms spectra 
always show an “inverted” shape (as typically seen in the RXTE 
data, e.g. Belloni et al. 2011), with variability power decreasing as 
a function of energy. We conclude that a relatively high fraction 
of disc flux is variable during the hard state of GX 339-4. How¬ 
ever, the increasing fraction of constant disc flux with luminosity 
in the soft band is responsible for the net drop of fractional rms as 
compared to power law-dominated energy bands. We argue that this 
trend continues throughout the outburst, thus leading to the disc be¬ 
ing mostly constant during the soft state, and significantly different 
from the disc seen in hard states. 


5.2 The scaling of the high-frequency PSD with luminosity in 
BHXRBs and AGN 

A number of correlations link the high-frequency variability prop¬ 
erties of BHXRBs and AGN, suggesting that the dominant vari¬ 
ability process at these frequencies is the same in both the classes 
of objects. In particular, considering a sample of 10 AGN and 2 
BHXRBs in the soft state, McHardy et al. (2006) found a linear 
scaling of the high-frequency break (a characteristic frequency of 
AGN PSDs, Edelson & Nandra 1999, commonly associated with Vh 
in BHXRBs) with the BH mass (Mbh) and the mass-accretion rate 
(liiEdd, i.e. normalized by the mass-accretion rate required to reach 
the Eddington luminosity), which defines the so-called ‘variability 
plane’. Koerding et al. (2007) found a similar scaling relation also 
for sources in the hard state (but offset with respect to the corre¬ 
sponding relation valid for soft-state sources). 


As a consequence, we would expect to observe v;, scaling 
linearly with the mass-accretion rate also in the hard-state obser¬ 
vations of GX 339-4 analysed in this paper. However, we see v;, 
scaling with LILgdd (with L = Ls-sorcv), with an average slope 
of ~ 0.21 ± 0.11, therefore not consistent with the results by 
McHardy et al. (2006) and Koerding et al. (20070 We note, how¬ 
ever, that our result is in agreement with the recent estimate ob¬ 
tained from the PSD analysis of a sample of 104 nearby AGN ob¬ 
served with XMM-Newton (Gonzalez-Martin & Vaughan, 2012), 
which reports a best-fit slope of 0.24 ± 0.28. This therefore sug¬ 
gests that, indeed, the same engine is at work in both AGN and 
BHXRBs, but that the dependence of the high-frequency variability 
on luminosity/mass-accretion is significantly shallower than previ¬ 
ously thought (McHardy et al. 2006; Koerding et al. 2007). We also 
investigated the dependence of the normalization of the PSD at the 
break, W,,, on the luminosity. We find that Ny^ decreases as a func¬ 
tion of LILEdd- Interestingly, this trend can counterbalance the v/, 
vs. LjLEdd scaling, producing the same high-frequency variability 
power at different luminosities. 

Previous studies have pointed out that the high-frequency PSD 
does not show a clear scaling with the luminosity/mass-accretion 
rate both in BHXRBs (Gierlihski et al. 2008) and AGN (O’Neil et 
al. 2005; Zhou et al. 2010; Ponti et al. 2012). We observe that the 
high-frequency fractional rms of GX 339-4 is constant in the hard 
and very hard band (dominated by the hard X-ray Comptonization 
continuum) at any luminosity probed. In particular we find that this 
quantity has a small scatter (4-10 percent, in agreement with re¬ 
sults by Munoz-Darias et al. 2011 and by Heil, Vaughan, & Uttley 
2012) in these two energy bands. The remarkable constancy of the 
high-frequency PSD provides us with a tool to measure the BH 
mass of accreting objects (Nikolajuk et al. 2004; Gierlinski et al. 
2008; Ponti et al. 2012; Kelly et al. 2011; La Franca et al. 2014). 
Indeed the high-frequency fractional rms (or equivalently the ‘ex¬ 
cess variance’, e.g. Nikolajuk et al. 2004; Ponti et al. 2012), while 
being independent of the mass-accretion rate, scales inversely with 
the BH mass, with a very small scatter (which mostly depends on 
the uncertainty on Mbh), of the order of less than 0.2-0.4 dex (e.g. 
O’Neill et al. 2005; Zhou et al. 2010; Ponti et al. 2012; Kelly et al. 
2011 ). 

As observed in Sect. 14.31 the dependences of v;, and are 
significantly different in the soft band. In particular, the constancy 
of the high-frequency fractional rms as a function of luminosity 
breaks down at low energies, with a scatter around the best-fit con¬ 
stant model increasing to ~ 30 percent. We ascribe this to the pres¬ 
ence of the additional contribution from the disc component. Inter¬ 
estingly the same behaviour is observed in AGN. Indeed, it is ob¬ 
served that the scatter in the relation between the short-time scale 
variability and the BH mass is larger in the optical band (domi¬ 
nated by the disc emission) than in the X-ray band (dominated by 
the power law), and the optical variability of AGN is significantly 
anticorrelated with the optical luminosity (Kelly et al. 2013). 

5.3 A possible interpretation for the offset of the 
high-frequency Lorentzian in the soft band 

The characteristic frequencies in the PSD of accreting systems 
(i.e. the peak-frequency of each Lorentzian) are likely related to 

^ Here we assume LboI L and LjjoilLEdd ^ hJ^Edd’ ^ “ 1 ~ 2, the 
exact value depending on the accretion efficiency of the flow (Fender & 
Munoz-Darias 2015). 
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fundamental time scales of the accretion flow at specific radii 
(e.g. Psaltis, Belloni & van der Klis 1999; Churazov, Gilfanov, 
& Revnivtsev 2001; Motta et al. 2014). Though not yet conclu¬ 
sive, the general picture that emerged over the years is that the 
low-frequency Lorentzian in the hard state of BHXRBs is asso¬ 
ciated with a characteristic time scale (most probably the viscous 
time scale) at the inner edge of the accretion disc (e.g. Churazov 
et al. 2001; Done et al. 2007). According to recent studies, varia¬ 
tions intrinsic to the seed-photons flux produce the long-time scale 
variability (i.e. the low-frequency Lorentzian) observed in the soft 
energy band (Sect. 14.1.11 Wilkinson & Uttley 2009; Uttley et al. 
2014). These variations are due to perturbations in the accretion 
rate, which travel inward and modulate the hard X-ray emission 
from the inner corona, thus producing the long-time scale variabil¬ 
ity observed also in harder energy bands (e.g. Lyubarskii 1997; 
Kotov, Churazov, & Gilfanov 2001; Arevalo & Uttley 2006; In¬ 
gram & van der Klis 2013). The existence of additional variability 
components at frequencies higher than those corresponding to the 
inner edge of the disc, can be explained provided the corona ex¬ 
tends down to the innermost stable circular orbit (ISCO, Motta et 
al. 2014) and is intrinsically variable on short time scales (e.g. Chu¬ 
razov, Gilfanov, & Revnivtsev 2001). This inner hot flow/corona is 
thus responsible for producing the high-frequency Lorentzian ob¬ 
served in hard and very hard energy bands. 

Such a high-frequency component is observed also in the soft band, 
and is linearly correlated with that in harder energy bands (as seen 
from covariance spectra, light grey triangles in Fig.|^ see also Ut¬ 
tley et al. 2011; De Marco et al. 2015). We showed that the high- 
frequency Lorentzian is significantly offset, by a factor of ~2, to¬ 
wards lower frequencies in the soft band with respect to higher- 
energy bands. The covariance spectra (light grey triangles in Fig. 
nil show that the high-frequency PSD contains significant contri¬ 
bution from variable disc emission, so the olfset, soft band high- 
frequency Lorentzian is likely linked to variability coming from the 
disc component. However, within the scenario depicted above, in¬ 
trinsic disc variability due to accretion rate perturbations occurring 
on time scales shorter than the viscous time scale at the inner edge 
of the disc is expected to be highly suppressed. Thus, a possible 
explanation is that this high-frequency variability is not intrinsic 
to the disc, but results from reprocessing of the short-time scale 
variable corona emission (e.g. thermal reprocessing; Wilkinson & 
Uttley 2009; Uttley et al. 2011; De Marco et al. 2015). In standard 
disc-corona geometries, the hard X-rays would irradiate a relatively 
extended region of the disc. Thus, the reprocessing region would act 
as a low-pass filter, damping the high-frequency variability power 
in the soft band (e.g. Arevalo & Uttley 2006; Uttley et al. 2014). 
As a consequence the soft band high-frequency Lorentzian would 
be narrower, and its peak-frequency would result slightly shifted 
towards low frequencies with respect to the corresponding peak in 
harder energy bands. This resembles the behaviour observed in the 
data of GX 339-4. According to this interpretation, the relative off¬ 
set of the peak-frequency of the high-frequency Lorentzian in the 
different bands should trace the variations of the geometry of the 
accretion disc/corona during the outburst. 
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